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Whether the pairing mechanism for superconductivity in iron-based superconductors is of itinerant 
or local character is a question still debated on. In order to investigate the influence of Fermi surface 
nesting, we calculate from first-principles the static susceptibihty of FeSe under applied pressures, 
both hydrostatic and non-hydrostatic. We show that the electronic band structures are highly 
sensitive to the way pressure is applied and confront our theoretical results with conclusions drawn 
from experiments. Given that the critical temperature of FeSe is quite universal as function of 
pressure, this is clear evidence that in FeSe the evolution of Fermi surface nesting cannot account 
for the evolution of the critical temperature with pressure. Hence we argue that the pairing in 
iron-chalcogenide compounds should be of strong-coupling spin fluctuation origin. 

PACS numbers: 74.70.Xa, 74.62. Fj, 74.20.Pq 



The discovery of a new class of iron based high-Tc su- 
perconductors in 2008 [T] has kicked off an avalanche of 
research in solid state physics, a lot of it dedicated to 
the key question of the superconducting pairing mecha- 
nism. A feature that has been agreed upon is that the 
parent compounds of the FeAs superconductors exhibit 
an antiferromagnetic (AFM) spin density wave (SDW) 
phase, as seen in LnOFeAs {Ln = La,Ce,Pr,Sm [2H1]) and 
^eFe2As2 (^e=Ba,Sr,Ca [MZl), or are on the verge of it, 
as in ^FeAs [A = Li,Na [8j|^). In such materials, su- 
perconductivity (SC) arises or is enhanced, respectively, 
when this long-range magnetic order is suppressed by 
doping or by application of physical or chemical pressure. 
In this respect, pnictide and chalcogenide compounds are 
much more flexible as compared to cuprate high-Tc su- 
perconductors, where the superconducting state can be 
reached only by doping, and intercalating layers between 
the copper-oxygen planes are therefore essential. This 
large flexibility has been underpinned by the discovery 
of SC in FeSe [T0Ul3) . a material lacking completely the 
intercalated planes of the iron-pnictides. This structural 
simplicity and its enormous response to external pressure 
(Tc increases from 8K at ambient pressure [T^ to 37 K 
at about 7 GPa made it a popular testing ground to 
study mechanisms of SC in iron-based superconductors. 

Theoretical studies have shown [TSHTT] that iron 
chalcogenides are signiflcantly stronger correlated than 
the iron pnictide materials. Since this goes hand-in-hand 
with the enhancement of localized magnetic moments, 
the applicability of the Fermi surface nesting (FN) sce- 
nario in this class of compounds has to be clarifled. In 



that scenario, SC can be regarded as a consequence of 
perturbations of the AFM ordering of the parent mate- 
rial [THH23] : and suppression of peaks in the susceptibility 
due to FN is an indicator for the occurrence of SC. 

In order to shed light on this issue we perform first- 
principle calculations of the static susceptibility of FeSe 
under pressure. Many independent experimental studies 
on pressure vs. Tc phase diagrams have been published 
for FeSe [TH [201 I^^H^^ . The general consensus seems to 
be that Tc starts around 8 K for ambient pressure, has a 
maximum at about 7 GPa and decreases again for higher 
pressures. While all the works mentioned above applied 
pressure hydrostatically, a recent study investigated how 
Tc changes under won-hydrostatical pressure [27j . It is 
interesting to observe that the corresponding phase dia- 
gram is quite similar, i.e. the way how pressure is applied 
seems to be of minor importance for SC. 

Our results show that suppression of FN cannot be 
reconciled with the occurrence of SC. Furthermore, the 
way how pressure is applied to the system substantially 
changes the electronic band structure and consequently 
the static susceptibility. Given the universality of Tc vs. 
pressure, this indicates that SC in FeSe is driven by short- 
ranged spin fluctuations between rather localized mag- 
netic moments [ ^ [ ^ 1^. 

Method: A very useful quantity for theoretical inves- 
tigations of the strength and changes of the FN is the 
static susceptibility matrix x'^(q) (in short x"). For a 
system of Bloch electrons - taking into account only the 
diagonal elements of this matrix - it reads in the random 
phase approximation 
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FIG. 1: Static susceptibility x° in the (Qx, Qi,, 0)-plane of FeSe 
at a hydrostatic pressure of 2.8 GPa. The vertical solid black 
lines mark the positions of the X points in the reciprocal unit 
cell and q^, qy are given in units 27r/a and 27r/&, respectively. 



where q represents a vector of the extended wave number 
space, n and m denote electron band indices, k belongs to 
the first Brillouin zone (BZ), and e„i(k) and /m(k) mean 
the energy dispersion of the mth band and the Fermi 
function, respectively. The reciprocal lattice vector G is 
defined such that k + q + G G BZ. We want to stress 
that the matrix elements (m, k|e~"'''|n, k + q + G) are 
fully included in our calculation. 

FN leads to a pronounced peak of x'^ at the nesting vec- 
tor q — qAT due to singularities of the integrand. Many 
undoped parent materials of iron-based superconductors 
show a marked FN and consequently a maximum of x'^ 
at qAT. A perturbation of the crystal structure due to ap- 
phed pressure or of the electronic band structure due to 
(electron or hole) doping is able to suppress the FN and 
the strong response of It is argued that the 

reduction of the collectively AFM-ordered electron spins 
in favor of the appearance of spin-flip processes may act 
as a generator for electron pairing and SC. 

All results presented in this work have been obtained 
on the basis of non-spinpolarized DFT FPLAPW band 
structure calculations using the WIEN2k package [31]. 
The exchange-correlation potential has been approxi- 
mated by the generalized gradient approximation GGA, 
and the base-centered orthorhombic crystal structure 
Cmma has been used. The numerical integration over 
the first Brillouin zone has been performed with a semi- 
analytical tetrahedron method. Each electron wave func- 
tion has been represented by about 4000 plane waves, 
and 52 electron bands, corresponding to an energy in- 
terval from — 1.2Ry below to 3Ry above the Fermi en- 
ergy, have been included. The momentum vectors q are 



FIG. 2: Static susceptibility in the {qx,0,0) direction of 
FeSe for different hydrostatic pressures {q^ is given in units 
27r/a). Apart from the vertical shift, due to the increasing 
density of states at the Fermi energy, x° shows very little 
dependence on pressure. 



expressed in units of (27r/a, 27r/6, 27r/c) throughout this 
work, where a, b and c are the FeSe unit cell lattice pa- 
rameters, and x° is presented in units (eV-ao)~^ with oq 
being the Bohr radius. The lattice and unit cell param- 
eters are taken from experiment, where we used data by 
Margadonna et al. [T3] and Kumar et al. [5T] for hydro- 
static pressures, and for the non-hydrostatic case data 
from Uhoya et al. [57] . 

Results: We start the discussion of the pressure depen- 
dence with the hydrostatic case. We checked that the 
results based on the crystal structures reported by Ku- 
mar et al. |21| are essentially the same compared to the 
results based on structures reported by Margadonna et 
al. [13] . In the following, we will hence show only results 
based on data by Kumar et al.. In Fig. [ija surface plot 
of the static susceptibility over the (g^;, gj,, 0)-plane for 
FeSe at a hydrostatic pressure of 2.8 GPa is presented. 
The dominating features are the marked peaks, indicat- 
ing a strong FN of the hole pockets at the F point and 
the electron pockets at the X points, corresponding to 
the FN vector qAr = (1,0,0). For a better quantita- 
tive discussion of the pressure dependence, we show a 
(ga;,0, 0) cut of the susceptibility in Fig. [2j In this dia- 
gram we compare the susceptibility in the (1, 0, 0) direc- 
tion for three structures at pressures 2.8 GPa (red, x), 
7 GPa (cyan, o) and 21 GPa (blue, A). Interestingly, 
the three calculated profiles have a quite similar shape. 
They differ only by a vertical shift, reflecting the increase 
of the electronic density of states N{e) at the Fermi en- 
ergy ep with increasing pressure, according to the rela- 
tion x°(0) = N{eF)/^Q where fio denotes the volume of 
the unit cell of the crystal. However, the prominent peak 
at (\N — (1,0,0) does not vanish with increasing pres- 
sure, leading to the conclusion that the strength of the 
FN is not pressure-dependent in the hydrostatic case. 
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FIG. 3: Static susceptibility x'^ in the (qx, 9;/, 0)-plane for 
FeSe at a non-hydrostatic pressure of 4 GPa. The vertical 
solid black lines mark the positions of the X points in the 
reciprocal unit cell and q^, Qy are given in units 27r/a and 
27r/6, respectively. 




FIG. 4: Static susceptibility x" in the (g^, gj,, 0)-plane for 
FeSe at a non-hydrostatic pressure of 12.9 GPa. The vertical 
solid black lines mark the positions of the X points in the 
reciprocal unit cell and q^, Qy are given in units 27r/a and 
27r/ti, respectively. 



The situation is different in the non-hydrostatic case, 
as demonstrated in Figs. [3] and |4j Here we present 
the static susceptibility for the compounds reported by 
Uhoya et al. H?] . Fig. [3] depicts the situation for 4 GPa, 
where the marked maxima at the X points, similar to 
the hydrostatic case, can still be observed. On the other 



FIG. 5: Static susceptibility in the {cix,Q,Q) direction for 
FeSe for different non-hydrostatic pressures (q^ is given in 
units 27r/a). shows a marked peak at the FN vector qjv for 
the 4 GPa case, which is still present, albeit reduced, for the 
7 GPa compound. At high pressures of 12.9 GPa no prominent 
feature can be seen. 



hand, for the high pressure case (12.9 GPa) shown in 
Fig. [4j these maxima are reduced significantly. This 
change of the susceptibility is even better demonstrated 
in Fig. [sj where we depict again a (ga;,0, 0) cut of the 
susceptibility. Here we compare for the case of non- 
hydrostatic pressures at 4 GPa (red, x), 7 GPa (cyan,o) 
and 12.9 GPa (blue, A). For the intermediate pressure 
regime of 7 GPa we interpolated the structural data from 
Uhoya et al. ^ for 4 GPa and 12.9 GPa. It is known 
from experiments [131 [51] [3S] that - in such moderate 
pressure regimes - the unit cell parameters change lin- 
early with pressure, supporting this approach. While the 
strong response at q^r is still present for the cases with 
4 and 7 GPa (albeit slightly reduced for the 7 GPa struc- 
ture), it has vanished completely for the high pressure 
case. This behavior of allows the conclusion that the 
FN in FeSe is rather strong for low and reduced consid- 
erably for higher non-hydrostatic pressures. 

In order to get better insight into this different FN 
behavior, we present in Figs. [6] and [7] the cuts of the 
Fermi surfaces in the {qx,qy, 0)-plane for low (4 GPa) and 
high (12.9 GPa) non-hydrostatic pressures. The different 
colors in these diagrams correspond to different bands, 
while the black dotted lines mark the edges of the unit 
cell in reciprocal space. 

For low pressures (Fig. [6]) , it can clearly be seen that 
there is a marked FN, originating mainly from the in- 
terband transitions between the hole-like pockets around 
the r point and the electron-like pockets around the X 
point. Different to doping, where hole pockets shrink 
and electron pockets grow - or vice versa - the effect 
of non-hydrostatic pressures is a bit more subtle: Large 
pressures lead to considerable changes in the a/b lattice 
parameters, and as a consequence, not only the sizes of 



4 



FIG. 6: Fermi surface qz = cut for FeSe at a non-hydrostatic 
pressure of 4 GPa. The different colors correspond to differ- 
ent electron bands that cross the Fermi energy and the black 
dashed lines mark the boundaries of the reciprocal unit cell. 




more detailed and reliable information about the position 
of the SDW-vector and the strength of the FN. This is 
based on the fact that involves the integration over 
the whole Brillouin zone and the inclusion of the matrix 
elements, as can be seen in Eq. [T] 

Summary: The results presented above have impor- 
tant implications concerning the behavior of SC in FeSe. 
While for low pressures the static susceptibility of this 
material exhibits a marked FN, the higher pressure prop- 
erties depend crucially on whether pressure is applied hy- 
drostatically or not. According to our calculations there 
is no visible reduction of the FN for applied hydrostatic 
pressures, even as large as 21 GPa. On the other hand, 
we observe that the FN is reduced significantly if pressure 
is applied non-hydrostatically. Based on these theoreti- 
cal results alone one might argue that - in the framework 
of FN-mediated pairing - SC in FeSe is only supported 
efficiently by non-hydrostatic pressures. However, such a 
result contrasts sharply with experiments, which report 
a considerable increase of Tc with increasing pressure in- 
dependent of the way pressure is applied. 

This leads to the conclusion, that SC in the FeSe com- 
pounds is not based on a weak-coupling nesting between 
electron and hole pockets, and therefore not mediated by 
spin fluctuations coming from an itinerant SDW state. 
Moreover, recent studies have shown [35H3E! that the pro- 
nounced increase of Tc with pressure cannot be explained 
by phononic Cooper-pairing, too. Therefore, staying in 
the paradigm of spin-mediated pairing, we would argue 
that the SC in FeSe is based on fluctuations between more 
localized magnetic moments. 



FIG. 7; Fermi surface = cut for FeSe at a non-hydrostatic 
pressure of 12.9 GPa. The different colors correspond to dif- 
ferent electron bands that cross the Fermi energy and the 
black dashed lines mark the boundaries of the reciprocal unit 
cell. 



the pockets, but also their shapes are altered substan- 
tially (Fig. [7|. In the case of low pressures the elec- 
tron and hole pockets are in good approximation cir- 
cularly shaped. For high pressures on the other hand, 
the pockets at the X points disappear completely and 
the Fermi surfaces at the F point are distorted to more 
complicated elliptical shapes. This results in a complete 
loss of the FN vector qjv and further supports our con- 
clusion drawn from Figs. |3][5] namely, that an increas- 
ing non-hydrostatic pressure causes a signiflcant decrease 
of the FN: While in the 4 GPa case the nesting vector 
Qtv = (liO,0) can clearly be observed, such a feature is 
not visible for a pressure of 12.9 GPa. We want to stress 
that the calculation of the complete susceptibility offers 
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